Cosmology with positive and negative exponential potentials 
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We present a phase-plane analysis of cosmologies containing a scalar field <f> with an exponential 
potential V oc exp( — \K(f>) where ft 2 = 8tyG and V may be positive or negative. We show that power- 
law kinetic-potential scaling solutions only exist for sufficiently flat (A 2 < 6) positive potentials or 
steep (A 2 > 6) negative potentials. The latter correspond to a class of ever-expanding cosmologies 
with negative potential. However we show that these expanding solutions with a negative potential 
are to unstable in the presence of ordinary matter, spatial curvature or anisotropic shear, and 
generic solutions always recollapse to a singularity. Power-law kinetic-potential scaling solutions are 
the late-time attractor in a collapsing universe for steep negative potentials (the ekpyrotic scenario) 
and stable against matter, curvature or shear perturbations. Otherwise kinetic-dominated solutions 
are the attractor during collapse (the pre big bang scenario) and are only marginally stable with 
respect to anisotropic shear. 
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I. INTRODUCTION 

Scalar fields are ubiquitous in modern cosmology 
whether part of the matter sector, such as the Higgs field, 
or dilaton or moduli fields in the gravitational sector. 
They extend the range of qualitative behaviour found 
in homogeneous and isotropic Friedmann-Robertson- 
Walker (FRW) models compared with conventional per- 
fect fluid models. Scalar fields with positive potential 
energy can drive an accelerated expansion and thus pro- 
vide possible models for inflation in the early universe or 
dark energy at the present epoch. 

In this paper we will consider models where the scalar 
field, <f>, has a simple exponential potential 



V = Vq cxp(—\K(, 



(1) 



where k 2 = 8itGn is the gravitational coupling and A is 
a dimensionless constant characterising the slope of the 
potential. Such potentials arise naturally in supergravity 
models The possible cosmological role of exponen- 
tial potentials has been investigated before as a means of 
driving a period of early inflation j§-§h ekpyrotic or pre 
big bang collapse ||-|)) or quintessence po|-p^| . Although 
a period of accelerated expansion requires V > we will 
also consider the case where V < 0. Negative poten- 
tials have recently been considered both in an expanding 
universe Q and in a collapsing phase 

The scale-invariant form makes the exponential po- 
tential particularly simple to study analytically. There 
are well-known exact solutions corresponding to power- 
law solutions for the cosmological scale factor a oc t p 
in a spatially flat FRW model but more generally 
the coupled Einstein-Klein-Gordon equations for a sin- 
gle field can be reduced to a one-dimensional system 
which makes it particularly suited to a qualitative anal- 
ysis fflJH llMpf. Here we adopt a system of di- 



mensionless dynamical variables jl5(] which have previ- 
ously been used to study positive exponential potentials 
in isotropic jl6],[l8) and anisotropic [|l7| cosmologies, and 
mostly recently the brane- world scenario |is| ]. 

In Section II we review the known results for the one- 
dimensional system with V > and extend the analysis 
to include negative potentials. We then describe in sec- 
tion III the effect of introducing an additional term in the 
Friedmann equation to describe the effect of introducing 
a barotropic fluid, spatial curvature, or anisotropic shear 
in a Bianchi type I cosmology, thereby raising the system 
to two-dimensional phase-plane. We discuss the implica- 
tions of our results for the stability of particular solutions 
in section IV, epsecially those solutions invoked in the 
epyrotic and other pre big bang scenarios. 



II. ONE-DIMENSIONAL PHASE-SPACE 

We consider first a scalar field with an exponen- 
tial potential energy density evolving in a spatially- 
flat Friedmann-Robertson- Walker (FRW) universe. The 
Klein-Gordon equation for the scalar field in an FRW 
model is 



3H<j> + 4t = , 



(2) 



where the Hubble parameter, H, is determined by the 
Friedmann constraint 



H z = — [\$ 2 + V 



(3) 



A homogeneous scalar field has energy density — 
(p 2 /2 + V and pressure = 2 /2 — V. We introduce 
dimensionless variables 



1 



K(f> 



y 



<W\v\ 

V3H 



(4) 



in which the Friedmann constraint (|3|) takes the simple 
form 



2 i 2 

x ±y 



= 1 



(5) 



Throughout we will use upper/lower signs to denote the 



two distinct cases of ±V > 0. 



measures the contri- 



bution to the expansion due to the field's kinetic energy 
density, while ±j/ 2 represents the contribution of the po- 
tential energy. Their ratio determines the equation of 
state 



P. 



v x 2 =f y 2 



V % ± y 2 



(6) 



For a positive potential this which ranges between a stiff 
fluid with P^ = when the kinetic energy dominates, 
and Ptf, = —p<j>, when the potential dominates. For a 
negative potential we have an ultra-stiff fluid with P^ > 
P<t>- 

We will restrict our discussion of the existence and sta- 
bility of critical points to the upper half plane y > 0, i.e., 
expanding cosmologies with H > 0, but the trajectories 
are symmetrical under time reversal, H — > —H, corre- 
sponding to reflection symmetry x — > x and y — » —y. 
Thus early time solutions in an expanding universe are 
the same as to the late time limit in a collapsing cos- 
mology. Without loss of generality we will consider only 
A > as the system is also symmetric under A — > — A 
and x — ► —x. 

The evolution equation (^) can then be written as an 
autonomous system ]l5|,^6| : 



x' = -3x(l - x 2 ) ± AW -y 2 



y' = xy \3x - \\ 



(7) 

(8) 



where a prime denotes a derivative with respect to the 
logarithm of the scale factor, N = ln(a). The constraint 
(0) then reduces the system to a one-dimensional phase- 
space corresponding to the unit circle for V > or hy- 
perbola for V < 0. 



A. Critical points 

Critical points correspond to fixed points (xi, yi) where 
both x' = and y' = 0. These are self-similar so- 
lutions as dimensionless quantities are invariant under 
time-translation N — > N + AN. For example we have 
H/H 2 =constant, where 



H 

H 2 



= -3x z 



(9) 



Integrating this equation shows that all critical points, 
where x — Xi =non-zero constant, correspond to a power- 
law solution for the scale factor in terms of cosmic time: 



a cx \t\ p , where p 



1 

3x1 



(10) 



The system (0) and (||) has at most three fixed points: 

A± Two kinetic- dominated solutions exist for any form 
of the potential with y^ = and xa + = +1 or 
xa_ = — 1. These are equivalent to stiff-fluid dom- 
inated evolution with a cx t 1 / 3 irrespective of the 
nature of the potential. This is the form of cosmo- 
logical evolution usually invoked approaching the 
singularity in the Einstein conformal frame in the 
pre big bang scenario pFj] , 

B A potential-kinetic-scaling solution exists for ±(6 — 
A 2 ) > with 



x B 



71 



ys = V 1 - y 



(11) 



This solution only exists for sufficiently flat poten- 
tials (A 2 < 6) for positive potentials, or steep po- 
tentials (A 2 > 6) for negative potentials. 

The power-law exponent, p — 2/A 2 , depends on the 
slope of the potential. For positive potentials with 
A 2 < 2 and H > this solution corresponds to 
the well-known power-law inflation solutions with 
an accelerating scale factor a > §. For negative 
exponential potential with A ^> 1 and H < this 
yields the accelerated collapse with_£ <C 1 recently 
invoked in the ekpyrotic scenario 



B. Stability 

Linear perturbations x — * x + u about the points 
xa + = +1 and xa- — — 1 have exponential be- 
haviour u oc e mN where m+ = y/6(\/E — A) and 
m_ = \/6(^/6 + A) respectively. Thus for positive 
A, xa_ = — 1 is always unstable and xa + = +1 is 
stable for A 2 > 6 but unstable for A 2 < 6. 



B Linear perturbations about the potential-kinetic 
scaling solution decay with eigenvalue m — (A 2 — 
6)/2 and this solution is thus always stable when 
this point exists for a positive potential, but unsta- 
ble when this solution exists for a negative poten- 
tial. 



C. Behaviour at infinity 

For positive potentials the one-dimensional system is 
compact, but for V < the trajectories can go out to 
infinity with x/y — ► 1 and the system reduces to 
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FIG. 1. One-dimensional phase-space for flat positive po- 
tentials, A 2 < 6. Arrows indicate evolution in cosmic time, 
t. Note that in the lower half-plane, H < 0, this has the 
opposite sense to A^ = ln(a). 

x' -> 3x 3 (12) 

so trajectories always approach infinity with x 2 = 
1/6(AT» — N). Equation (||) can then be integrated to 
obtain a solution for H and hence the scale factor in 
terms of proper time 

A/*-7Voc {t-U) 2 , (13) 

showing that the behaviour at infinity always corresponds 
to recollapse, i.e., a maximum of the expansion at finite 
time t*. Thus expanding solutions in the upper half- 
plane are linked to collapsing solutions in the lower half- 
plane at infinity. When H changes sign, both x and y 
also change sign, so solutions in Figures || and [| exiting 
top-right are linked to solutions bottom-left and similarly 
top- left to bottom-right. 

D. Qualitative evolution 

a. Flat positive potentials (V > 0, A 2 < 6) . All 
three critical points exist. Point B is the stable late- 
time attractor and points A— and A + are unstable repel- 
lors. Hence generic solutions start in a kinetic-dominated 
regime and approach the kinetic-potential scaling solu- 
tion at late times. 

b. Steep positive potentials (V > 0, A 2 > 6) . Only 
the two kinetic-dominated fixed points A + and A- exist. 
For A > 0, generic solutions start at xa_ = —1 and 
approach xa + = +1 at late times. 

c. Flat negative potentials (V < 0, A 2 < 6) . Only 
the two kinetic-dominated fixed points A + and A_ exist 
and both are unstable repellors. Thus generic solutions 
begin in one of the kinetic dominated solutions and go 
out to infinity, i.e., recollapse. Thus these cosmological 
solutions have a finite lifetime, beginning in a kinetic- 
dominated big bang (point A- or A + ), but recollapsing 
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FIG. 2. One-dimensional phase-space for steep positive 
potentials, A 2 > 6. 




FIG. 3. One-dimensional phase-space for flat negative po- 
tentials, A 2 < 6. 

after finite time and collapsing to a big crunch (at A + or 

d. Steep negative potentials (V < 0, A 2 > 6) . All 
three critical points exist. For A > 0, point A + is 
the stable late-time attractor. Thus left hyberbola de- 
scribes solutions that leave xa_ = — 1 and go to infin- 
ity. Generic solutions on the right hyperbola start in the 
kinetic-potential scaling solution at point B and then ci- 
ther become kinetic-dominated (approaching A + ) or go 
to infinity (recollapse). 

III. TWO-DIMENSIONAL PHASE-SPACE 

We now consider the effect of an additional component, 
not explicitly coupled to the scalar field, but gravitation- 
ally coupled via the Friedmann constraint for the Hubble 
expansion: 

H2 = Y (^ 2 + v + p^ ■ ( 14 ) 
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Potential 


Evolution 


Positive 


flat 


A- 




B+ 






big bang 


expansion 


future infinity 






A+ 


_> 


B+ 






big bang 


expansion 


future infinity 






B- 




A- 






past infinity 


collapse 


big crunch 






B- 


-» 


A+ 






past infinity 


collapse 


big crunch 




steep 


A- 


-> 


A+ 






big bang 


expansion 


future infinity 






A+ 


-> 


A- 






past infinity 


collapse 


big crunch 


Negative 


flat 


A- 




A+ 






big bang 

o o 


expansion and recollapse 


big crunch 






^+ 




A- 






big bang 


expansion and recollapse 


big crunch 




steep 






A+ 






big bang 


expansion 


future infinity 










A- 






big bang 


expansion and recollapse 


big crunch 






A_ 




B- 






big bang 


expansion and recollapse 


big crunch 






A+ 




B- 






past infinity 


collapse 


big crunch 



TABLE I. The evolution of one-dimensional systems. 




FIG. 4. One-dimensional phase-space for steep negative 
potentials, A 2 > 6. 



where p 1 obeys the continuity equation 

Pi 



-3# (p 7 + P 7 ) 



(15) 



and the scalar field obeys the same Klein-Gordon equa- 
tion (|). 

We will consider the simplest case of a barotropic fluid, 
with pressure P 7 = (7 — l)/) 7 where 7 =constant and 
< 7 < 2 for conventional fluids (such as dust, 7=1, 



or radiation, 7 = 4/3). In this case we again have an 
autonomous system pi| in terms of dimensionless vari- 
able x and y defined in Eq. (Q) and a new dimensionless 
variable w 



(16) 



yp 7 

V3H ' 



The evolution equations (|^) and ( |T 5 |) can then be writ- 
ten as an autonomous system |l5| , |l6| |: 



3x 1 - x 1 



1 U 2 , 3 7 2 
/ = y I 6x + —w 



(-7 + 2x 2 + jw 2 ) 



(17) 
(18) 
(19) 



where as before a prime denotes a derivative with respect 
to N. 

The Friedmann constraint (Q) takes the form 



2,2, 2 

2; ± y + w 



1 



(20) 



and reduces the system to a two-dimensional phase-space 
corresponding to the unit sphere for V > or hyper- 
boloid for V < 0. Note that w = remains an invariant 



4 



one-dimensional subspace. We will restrict our analysis 
to expanding cosmologies, H > 0, corresponding to the 
upper-quadrant y > and w > 0, noting that the sys- 
tem is symmetric under t — > —t, H — > —H, y — > — y and 
w — * —w. 

The same system of equations can also be used to de- 
scribe Bianchi type I cosmologies |oL where the shear 
appears in the constraint equation (jlj) like a fluid with 
7 = 2, or spatially curved FRW models (lq , where the 
curvature appears in the Friedmann equation like a fluid 
density with 7 = 2/3 (and p 1 < 0, for positive curved 
space, or p 7 > 0, for negative curvature). In particu- 
lar the stability analysis for the barotropic fluids with 
7 = 2 or 7 = 2/3 is also applicable to scalar fields with 
exponential potentials in Bianchi type I or curved FRW 
models, where the flat FRW models are an invariant one- 
dimensional subspace. 



A. Critical points 

Critical points correspond to fixed points (xi, yi, tUj) 
where x' = 0, y' = and w' = 0, and these are self- 
similar solutions with, for example, 



H 



3 7 



(21) 



This corresponds to a power-law solution for the scale 
factor 



a oc \t\ p , where p 



6xf + 3jwf 



(22) 



The system (pf), (ph and (|) has at most five fixed 
points given by the vacuum/scalar field solutions A±, B, 
and two new points: 



C A fluid- dominated solution 

xc = , yc = , u>c = 1 , 



(23) 



always exists, for any form of potential, correspond- 
ing to a power-law solution given by Eq. (J22j) and 
p = 2/3 7 . 

Z) A fluid-potential-kinetic-scaling solution exists for 
steep positive potentials: 



and is independent of the slope of the potential. 
This is the classic example of 'tracking' behaviour 
which occurs for a wide variety of steep positive 
potentials 



B. Stability 

The stability of A± and B with respect to perturba- 
tions in the w = sub-space is unchanged, yielding one 
eigenmode and eigenvalue as discussed in sect. [II B| . In- 
troducing fluid perturbations yields a second eigenmode 
w' = 3{2x\ ! -7)«;/2: 

A± Fluid perturbations evolve as w oc e mN where the 
eigenvalue m = 3(2 — 7)/2. Thus the kinetic- 
dominated solutions are always unstable to fluid 
perturbations for 7 < 2. For a stiff fluid with 7 = 2 
(or Bianchi type I shear) the kinetic-dominated so- 
lutions are marginally stable. 

B Fluid perturbations evolve as w oc e mN where the 
eigenvalue m = (A 2 — 37)/2. Thus the kinetic- 
potential scaling solution is stable with respect to 
fluid perturbations for sufficiently flat potentials, 
A 2 < 37, but unstable for A 2 > 37. In particular, 
the kinetic-potential scaling solutions in spatially- 
flat FRW models are stable with respect to spatial 
curvature for A 2 < 2, yielding with power-law in- 
flation for positive potentials W. 

The two new points also have two eigenmodes: 

C The fluid-dominated solution has two eigenmodes 
x' = -3(2 - j)x/2 and y' = fryy/2. Thus the so- 
lution is stable to kinetic energy perturbations for 
7 < 2, but unstable to potential energy perturba- 
tions for 7 > 0. Thus for < 7 < 2 it is a saddle 
point and hence unstable to generic perturbations. 

D Because point D can only exist for positive po- 
tential the phase-space reduces to that studied in 
Refs. |l6| where it was shown that this scaling so- 
lution is always the late-time stable attractor when 
it exists. 



C. Behaviour at infinity 



X D 



*1 

2A ' 



3 (2-7)7 

37 



w D = l-^. (24) 



This solution only exists for potentials that are 
both positive (V > 0) and sufficiently steep (A 2 > 
37 2 ). The power-law exponent in Eq. (B2n, p = 
2/37, is identical to that of the fluid-dominated so- 
lution and depends only on the barotropic index 



For positive potentials the one-dimensional system is 
compact, but for V < the trajectories can go out to 
infinity. We have 



= -3- 



y 



(25) 



and hence as y — > 00 we find that generic solutions ap- 
proach x/y — > — 1. 

The remaining equations can then be written as 



5 




FIG. 5. Two-dimensional phase-space for flat positive po- 
tentials, A 2 < 37. 

(* 2 + \u?y 3(x 2 + If (26) 

so trajectories approach infinity with x 2 + (7/2) w 2 — 
1/6 (iV* - N). Equation (0) can be integrated to obtain 
a solution for the scale factor in terms of proper time 

N*-N(x (t-U) 2 , (27) 

showing that the behaviour at infinity always corresponds 
to recollapse, i.e., a maximum of the expansion at finite 
time i* . Thus expanding solutions in the upper quadrant 
(w, y > 0) are linked to collapsing solutions in the lower 
quadrant (w,y < 0) at infinity. When H changes sign, 
x, y and w all change sign. 

D. Qualitative evolution 

e. Flat positive potentials (V > 0, X 2 < 3~f) . Four 
critical points exist in the upper quadrant y > 0. Points 
A + and A- are unstable repellors, C is the unstable fluid 
dominated solution, while B is the stable late-time at- 
tractor. Generic solutions begin at one or other of the ki- 
netic dominated solutions {A + or A J) and approach the 
potential-kinetic solution (B) at late times. See Figure ||. 
They can approach arbitrarily close to the fluid domi- 
nated saddle point before reaching the potential-kinetic 
solution. 

/. Intermediate positive potentials (V > 0, 37 < A 2 < 
6) . All five critical points exist for y > 0. Generic so- 
lutions begin kinetic dominated ( A + or A_ ) and may ap- 
proach the fluid dominated (C) or potential-kinetic (B) 
saddle points before going to the fluid-potential-kinetic 
scaling solution (D) at late times. See Figure |(| 

g. Steep positive potentials (V > 0, X 2 > 6) . Four 
critical points exist. Points A + and A_ are unstable re- 
pellors, C is the unstable fluid dominated solution, while 
D is the stable late-time attractor. Solutions can be- 
gin in the kinetic dominated solutions and may approach 




FIG. 6. Two-dimensional phase-space for intermediate 
positive potentials, 37 < A 2 < 6. 




FIG. 7. Two-dimensional phase-space for steep positive 
potentials, A 2 > 6. 

the fluid dominated saddle point (C) before going to the 
fluid-potential-kinetic scaling solution (D) at late times. 
See Figure 0. Those that start out kinetic dominated can 
approach the scaling solution via the fluid dominated so- 
lution. 

h. Flat negative potentials (V < 0, X 2 < 6) . Three 
critical points exist. Points A + and A- are unstable re- 
pellors and C is the unstable fluid dominated solution. 
Generic solutions begin in one of the kinetic dominated 
solutions (A + or A J) and go to infinity along the line 
x = —y, before recollapsing towards the kinetic domi- 
nated collapse solution A + . See Figure ra. Expanding 
cosmologies that start out kinetic dominated can reach 
arbitrarily close to the fluid-dominated saddle point be- 
fore recollapsing. Generic solutions have finite lifetime 
before recollapsing to a singularity. 

i. Steep negative potentials (V < 0, A 2 > 6) . Four 
critical points exist and all are unstable. Generic solu- 
tions begin either kinetic dominated (A- for </> < 0) or 
kinetic-potential scaling (B for <fi > 0) and go to infinity 
along the line x = —y, before recollapsing towards the 
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FIG. 8. Two-dimensional phase-space for flat negative 
potentials, A 2 < 6. 




FIG. 9. Two-dimensional phase-space for steep negative 
potentials, A 2 > 6. 

kinetic-potential collapse solution (B) . See Figure [| If 
an expanding cosmology begins in the kinetic dominated 
solution (A_ ) it may approach the fluid dominated sad- 
dle point (C) before going to infinity. If an expanding 
cosmology begins in the potential-kinetic solution (£?) it 
may approach both the kinetic dominated saddle point 
(A + ) and the fluid dominated saddle point (C) before 
going to infinity. Generic solutions have finite lifetime 
before recollapsing to a singularity. 

IV. DISCUSSION 

We have presented a phase-space analysis of the qual- 
itative evolution of cosmological models with a scalar 
field, (p, with positive or negative exponential potential, 
V cx exp(— Xncj)). For positive potentials that are suffi- 
ciently flat (A 2 < 6) there is a self-similar scaling solution 
where potential and kinetic energies remain proportional. 
This is the basis of the well-known power-law inflation- 
ary solutions Q with a cx t p and p > 1/3. There is 
no such scaling solution for flat negative potentials, but 



sufficiently steep (A 2 > 6) negative potentials do have a 
self-similar power-law solution with p < 1/3. 

It has recently been found in the case of polynomial 
potentials |^2| that expanding cosmologies with a neg- 
ative potential always recollapse. By contrast, we find 
a expanding solutions with a sufficiently steep (A 2 > 6) 
negative potential that start close to the scaling solution 
at early times but approach the kinetic-dominated solu- 
tion at late times where the potential is negligible and 
expansion can continue indefinitely. Cosmologies with 
flat (A 2 < 6) negative potentials do always recollapse. 

The scaling solution is always the late-time attractor 
for flat positive potentials in an expanding universe, but 
the scaling solution with steep negative potentials is al- 
ways unstable in an expanding universe. Conversely the 
kinetic-dominated solution with stiff equation of state 
{P(j> — P<f>) is always the attractor approaching a big 
crunch singularity in a collapsing universe with positive 
potential, but for steep negative potentials the attrac- 
tor solution is the scaling solution with an ultra-stiff 
(P$ > p^) equation of state. This is quite contrary to 
the usual expectation that any scalar field potential be- 
comes negligible as the singularity is approached. For 
steep negative potentials the isotropic ultra-stiff scaling 
solution remains the attractor even in the presence of 
conventional matter perturbations, spatial curvature and 
anisotropic shear. 

It is interesting to apply these results to recently pro- 
posed cosmological models based on exponential poten- 
tials in a 'pre big bang' phase. The ekpyrotic model 
has a steep negative potential (A 2 ^> 6). Thus the 
kinetic-potential scaling solution is the attractor solu- 
tion approaching the big crunch singularity in a collaps- 
ing universe. It is stable to both conventional matter 
perturbations, curvature and shear perturbations. This 
echoes recent results found for anisotropic cosmologies 
in Randall-Sundrum-type brane- world scenarios p3| , but 
is in contrast to the pre big bang scenario for a mass- 
less dilaton field |20|| which is only marginally stable to 
anisotropic shear p4j. 

Another topical example is Finelli and Branden- 
berger's collapsing model || that gives a scale-invariant 
spectrum of comoving curvature perturbations during the 
collapse phase. This corresponds to a kinetic-potential 
scaling solution with p = 2/3, i.e., A 2 = 3, so requires 
a flat, positive potential. Thus this solution is unstable 
in a collapsing universe where generic solutions approach 
the kinetic-dominated regime. 
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